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We develop a theory for spin nematic ordering at finite temperatures in three-dimensional spatially 
anisotropic magnets consisting of weakly coupled frustrated spin--| chains. This theory is applied 
for weakly coupled J1-J2 chains with ferromagnetic nearest-neighbor J\ and antiferromagnetic next- 
nearest-neighbor J2 in magnetic fields. Combining the field theory technique with density-matrix 
renormalization group results, we complete the finite-temperature phase diagram in magnetic fields, 
which possesses spin bond-nematic and incommensurate spin-density-wave ordered phases. Effects 
of a four-spin coupling are also studied. The relevance of our result to quasi-one-dimensional edge- 
shared cuprate magnets such as LiCuV04 is discussed. 
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Introduction.— Multipolar orders have long been stud- 
ied in the context of heavy fermion systems. In recent 
years, their possible realization has been discussed much 
actively in quantum spin systems [l|-[3 • N-th order mul- 
tipolar order parameters are here defined by the expec- 
tation values of symmetrized tensor products of N spin 



operators S ri , S r 



In the N-th multipolar or- 



dered states, the 7V-th multipolar order parameter has a 
finite value, while all of the M(< 7V)-th multipolar order 
parameters vanish. Among them, a quadrupolar state is 
called a spin nematic state, which has a finite quadrupo- 
lar order (S+S^ + h.c.) 7^ and zero dipolar moment 
(S r ) = 0. In usual magnets, however, a spin order with 
a finite local magnetization (S r ) occurs at sufficiently 
low temperatures. Geometrical frustration is hence ex- 
pected to be an important ingredient for the emergence 
of multipolar order [l[. 

In the spin- \ systems, multipolar operators cannot be 
defined in a single site because of the commutation rela- 
tion of spin-^ operators. They reside on bonds between 
different sites [JEl], which is a significant difference from 
the multipolar phases in heavy fermion or higher-spin 
systems [7[. From this property, it is generally hard 
to develop controllable, effective theories for multipolar 
phases in spin-^ magnets compared to those in other sys- 
tems. Hence, most of the theoretical studies for multi- 
polar phases in two- and three-dimensional (3D) spin- 
| systems have been done by numerical computations 
for finite-size systems |l|, and analytic theories for the 
multipolar phases, especially, in wide temperature and 
external-field ranges, have been less developed. Mean- 
field theories have been developed quite recently for the 
spin nematic ground state [H, [s| . 

In ID spin-^ systems, on the other hand, various 
powerful theoretical/numerical techniques are applica- 
ble. Thanks to them, it has been shown recently 
that, in applied external magnetic field, three multipolar 



Tomonaga-Luttinger (TL) liquid phases 0, Q emerge in 
the spin-^ J1-J2 chain with ferromagnetic (FM) nearest- 
neighbor coupling Ji < and antiferromagnetic (AF) 
next-nearest-neighbor one J 2 > 0, whose Hamiltonian is 
given as 



(1) 



n=l,2 j 



Here Sj is the spin-^ operator on site j, and H is the 
external field. Near the saturation, quadrupolar (ne- 
matic) S^S^ +1 , octupolar S^ +1 S^ +2 i anc ^ hexadecap- 
olar Sj 1 S^ +1 S^ +2 Sj+3 operators exhibit quasi long-range 
order for the range - 
-2.7, and -3.76 < 



2.7 < Jx/Ja < 0, -3.5 < J ± / ' J 2 < 



J1/J2 



< 



-3.5, respectively, while 



the transverse spin correlator (£■ 5^) decays exponen- 
tially due to the formation of multiple-magnon bound 
states Q. These multipolar TL liquid phases expand 
down to a low- field regime, where the dominant correla- 
tion turns to an incommensurate longitudinal spin den- 
sity wave (SDW) type. 

The J1-J2 spin chain is expected to be an effective 
model for a series of quasi- ID edge-shared cuprate 
magnets such as LiCuVO^JTMl , Rb 2 Cu 2 MosOi 2 0, 
PbCuS0 4 (OH) 2 H3, E37 LiCuSb0 4 Q and 
LiCu 2 2 [20[. The experimentally estimated cou- 
pling ratio J1/J2 for LiCuV04 is well inside of the 
spin nematic TL-liquid phase in the Ji-J 2 chain. These 
theoretical and experimental results have motivated 
further searches of the spin nematic quasi and true 
long-range ordered phases in the high-field regime of 
LiCuV04 11, 15[. In addition, recent experiments in the 
intermediate field regime found incommensurate SDW 
oscillations 12-liJ whose wave vector nicely agrees 



with the theoretical prediction for the TL liquids of 
two-magnon bound states 0,11. It is still obscure how 
3D spin nematic and SDW ordered phases are induced 
with lowering temperature in quasi- ID magnets with 
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FIG. 1: (color online) Spatially anisotropic spin model con- 
sisting of weakly coupled spin-^ J1-J2 chains. We introduce 
inter-chain couplings J yi , J V2 , J V3 in the x-y plane. Similarly, 
J Zl , J Z2 5 Jz 3 are present in the x-z plane. 



weak interchain couplings, and how both of them are 
described in a unified way. 

In this paper, we develop a general theory for spin 
nematic and incommensurate SDW orders in spatially 
anisotropic 3D magnets consisting of weakly coupled J\- 
J2 spin chains in a wide magnetic-field range. Making 
use of field theoretical and numerical results for the J\- 
J2 spin chain, we obtain finite temperature phase di- 
agrams, which contain both spin nematic and incom- 
mensurate SDW phases at sufficiently low temperatures. 
From them, we reveal some characteristic features in the 
ordering of weakly coupled J\- J 2 chains. We also discuss 
the relevance of our results to real compounds such as 
LiCuV0 4 . 

Model.— Now, we start with the definition of our model 
of spatially anisotropic magnets shown in Fig. HJ The 
corresponding Hamiltonian is expressed as 



%3D — Hr + ^int? 



(2) 



where r = (r yi r z ) denotes the site index of the square 
lattice in the y-z plane, 7~L r denotes the Hamiltonian (pp) 
for the r-th J1-J2 chain along the x axis in magnetic 
field and Hint is the inter-chain interaction. In Hint, 
we introduce weak inter-chain Heisenberg-type exchange 
interactions with coupling constants J Vi and J z . (i = 
1, 2, 3) defined in the x-y and x-z planes, respectively (2lj . 

Spin-7^ J1-J2 chain.— Under the condition \ J yi , Zi \ *C 
I «7i,2 1 ? decoupled J1-J2 spin chains l~L r may be chosen 
as the starting point for analyzing the 3D model T~Lsd- 
Static and dynamic properties of the multipolar TL 1L 
uids in the J1-J2 chain H r have been well studied 
Q. The low-energy effective Hamiltonian for the spin 
nematic TL liquid phase (— 2.7 < J1/J2 < 0) is given by 



K ff = fdx t [KAd x o:) 2 + K-\d x r v ?] 



where x = a^j (the length ao of J\ bond is set equal to 
unity), (<fi±(x), 6±(x)) is the canonical pair of scalar bo- 
son fields, and v± and K± are, respectively, the excitation 
velocity and TL-liquid parameter of the (</>±,0±) sector. 
The term with the coupling G_ makes (/>_ pinned, induc- 
ing an excitation gap in the ((/>_,#_) sector. Physically, 
the gap corresponds to the magnon binding energy E^. 
On the other hand, the sector is described by 

a massless TL liquid. Vertex operators are renormalized 
as ( e ^V^+(*) e -^V^M0)} + = |2/ x |« 2 ^+/2 for | x | > 1? 

in which (• • • )± is the expectation value of the (4>±,0±) 
sector. 

Spin operators Sj jr are also bosonized as 

Sl r « M + d x (<l> r + + (-l^'O/V^ 

+ (-l) q A 1 cos[v^((/>; + + 27rMq] + • • • , 

(4a) 

5+ r «c < ^ F ^+(- 1 ^-){(-l)«B 

+ B x cos[v^(^; + + ZttMq] + •••}, (4b) 

where M = (Sj^), q = § (^~) f° r 3 = even (odd), 
and both A n and B n are nonuniversal constants. Utiliz- 
ing Eqs. (j3j) and (j4]), we can evaluate spin and nematic 
correlation functions at zero temperature (T = 0) as fol- 
lows 0,1,0: 



(5a) 

(5b) 
(5c) 



<S+So"> « Blco^j/2)(2/\j\fl^g_(x) 

{S Z 3 S Z Q ) « M 2 + (A 2 /2)|(e^-)_| 2 

xcos[^(M-l/2)](2/|j|)^/ 2 - 



-2/K+ 



+ G_ sin(TrM) sin( V4tt0!: + ttM) 



(3) 



where #_(x) = ( e W^M*) e W^-(o)}_ 5 c i s a constant 
and we have omitted the index r . The function g_ (x) de- 
cays in an exponential fashion like ~ x~ 1 / 2 e~ x ^~ . The 
TL-liquid parameter which is less than 2 in the low 
magnetization regime, monotonically increases as a func- 
tion of M [3[ and K+ 4 at the saturation. Thus, the 
spin nematic correlation is stronger (weaker) than the in- 
commensurate SDW correlation in high-field (low-field) 
regime with K+ > 2 (K+ < 2). 

The correlation length is related to v_ via v_ = 
(,-Eb under the assumption that the low-energy the- 
ory for the sector has Lorentz invariance. The 
velocity v+ has the relation v+ = 2K+/(irx), where 
X = dM/dH is the uniform susceptibility. The values 
of K+, and x are an determined with reasonable 
accuracy by using density-matrix renormalization group 
(DMRG) method 0, HH . Thus, v± can be quantitatively 
evaluated as depicted in Fig. O The figure shows that 
V- is always larger than and it is consistent with the 

perturbative formula v± « v(l ± KJi/fjrv) H ) in the 

weak IJ1I/J2 regime, in which u and if are respectively 
the spinon velocity and the TL-liquid parameter for the 
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FIG. 2: (color online) Magnon binding energy Eb (a)(c) and 
excitation velocities v± (b) (d) as a function of M in the spin 
nematic TL liquid phase in the spin-| J1-J2 chain at T = 0. 



single AF- J 2 chain. We also note that v+ approaches to 
zero in the vicinity of the saturation M — >• ~. 

Analysis of the 3D model— Let us now analyze the 
3D model ([2]) starting with the effective theory of the 
J1-J2 chain. We first bosonize all of the inter-chain 
couplings in Hint through Eq. (j4j). To obtain the 
low-energy effective theory for Eq. (|2]), we trace out 
the massive ((j) r _,6 7 l) sectors in the Euclidean action 
<?tot = <?o + <Smt Ym the cumulant expansion = 
So + (S int )- ~ \ ((5? t )_ - (5 int )2 ) + • • • , where <S and 
<Si n t are respectively the action for the TL-liquid part 
of the sectors and that for the inter-chain cou- 

plings. This expansion corresponds to the series expan- 
sion in J yiyZi /v-. The resultant effective Hamiltonian 
is expressed as H^q = Ho + Hsdw + %Ne + • • ■ • Here, 
no = ErI dxV ~Y [K + (d x r + ) 2 + K-\d x r+) 2 ] is the TL- 
liquid part, and Hsdw and ^Ne are, respectively, ob- 
tained from the first- and second-order cumulants as fol- 
lows: 



pairs of chains, where r' = r + e a (a = y,z), e a denotes 
the unit vector along the a- axis, and we have assumed 
that the field </> + smoothly varies in x. The first-order 
term Hsdw contains an inter-chain interaction between 
the operators e ±ly ^^+ , which essentially induces a 3D 
spin longitudinal order. Similarly, the term H^e contains 
an inter-chain interaction between the spin nematic op- 



erators sf r sf +lr 



(-iye ±W47v6 +, which enhances 3D 



spin nematic correlation. We should notice that the ef- 
fective theory W?Jl is reliable under the condition that 
temperature T is sufficiently smaller than the binding 
energy and the velocities v±. 

Both the couplings GsDW,Ne can be numerically eval- 
uated by using the correlation functions estimated with 
DMRG method 0, EH- ^sdw corresponds to the am- 
plitude of the leading term of the longitudinal correla- 
tor (SjSq) given in Eq. ([5]) and G^e can be evaluated as 

G Ne * ttvZ 1 Ef=iti/2) 1/K+ j(S?So) 2 - We have checked 
that the finite size-correction to the sum is small enough 
when the cut off L is larger than £_ . We emphasize that 
there is no free parameter in the field-theoretical Hamil- 
tonian H^ff. 

To obtain the finite-temperature phase diagram, we 
apply the inter-chain mean-field (ICMF) approxima- 
tion @, [IK to the effective Hamiltonian %l§ . To 
this end, we introduce the "effective" SDW operator 
C*sdw = e 27r ^2- M )i e ^^+ and the spin nematic oper- 
ator N e = (-l) j e iV ^ e +. Within the ICMF approach, 
the finite-temperature dynamical susceptibilities of O a 
(A =SDW or Ne) above 3D ordering temperatures are 
calculated as 



1 



4(% 



1D (k x ^Y 



where k = (k yi k z ) is the wave vector in the y-z plane, uo 
is the frequency, and the effective coupling constants 
are given by 

J e s ff DW (/t) =G S dw E [J^coska - J a2 sm(k a - ttM) 



, a3 sin(k a - 

Hsdw = Gsbw / y E E cos(v^(</>; - <)) J e N ff e (fc) =G Ne E [4 1 



• ttM)] , (8a) 
— (J a2 — Ja 3 ) 2 } cosk a . (8b) 



a = y,z 
(r f =r+e a ) 



a=y,z 



- j Q2 sin(V^(^; -</»;') -ttM) 

+ J a3 sin(V^(0; -</>;') + ttM)] , (6a) 

%Ne = G^e J — ^ [Jal 2 ~ (Ja2 ~ Ja3) 2 ] 



x cos(\/4^(6>; 



The ID susceptibilities x^°(fc a; ,w) = 

are analytically computed by using field theory tech- 
nique (P = l/T and e +0) [26]. Those for SDW and 
spin nematic operators respectively take the maximum 
value at /c™ ax = (\ - M)tt and tt; XsDw(^ ax ^°) = 



(6b) 



sm(^)2J(*± 1 



±^2 



and 



with coupling constants G S dw = Af|(e*^-) _| 2 Q 
and £?Ne = — 4^7- / dxv-drg-(x^r) 2 (r is imaginary 
time). The summations run over all nearest neighbor 



V-\- v / 

where B(x,y) is beta function. 

The transition temperature of each order is obtained 
from the divergent point of its own susceptibility at uo — >• 



0, which is given by 

l + Min fc [J e A ff (fe)] X ^(fcr x ,0) = 0. 



(9) 



The 3D ordered phase with the highest transition temper- 
ature T c is realized below T c . From this ICMF scheme, 
we can determine the phase diagram for with ar- 

bitrary combination of weak inter-chain couplings J yilZi - 
We should note that, when approaches to zero, the 
present framework becomes less reliable and we need to 
consider sub- leading terms in T~L^ '. 

From Eqs. (|5} and (|9]), we find that the ordering wave 
numbers k VjZ tend to be a commensurate value k VjZ = 
or 7r (see also the comment in Ref. I23I ). Thus the SDW 
ordered phase has the wave vector k x = {\ — M)tt and 
k VyZ = or 7r. This agrees with the experimental result 
in the intermediate-field phase of LiCuVCU [HI, [Hj . For 
the spin nematic ordered phase, we find the commensu- 



rate ordering vector (k x ,k y ( z )) 



(tt,0) for I J 



2/1(21) I 



> 



and Jy 3 are dominant, we find the similar tendency. The 
systems with dominant J V2 and J ys resemble the experi- 
mental proposal for LiCuVCU [1Q[ , where a new phase ex- 
pected to be a 3D nematic phase has been observed only 
near the saturation From the calculations for the 

cases of \J\\/ J2 — 0.5, 1.0, and 2.0, we find that the ne- 
matic phase region in the M-T phase diagram generally 
becomes smaller with increase in | J\ \ / J 2 since the value 
g~(x) in G^e decreases. The shrinkage of the nematic 
phase was also discussed in the ground state of multi-leg 
J1-J2 ladders 27]. When there is a certain frustration 
in interchain couplings, however, the nematic phase re- 
gion can expand, as shown in Fig. Efc). When the signs 
of J yi and J y2 (Jy 3 ) are opposite, the contribution to the 
coupling JffP w from the first-order perturbation becomes 
very weak, which expands the 3D nematic ordered phase 
down to a relatively lower-field regime. 

Effects of four-spin term.— Finally, we study effects of 
an inter-chain four-spin interaction on the phase diagram. 
The four-spin term we consider is 

Hi = _j 4 £ 5t p fi^ +liP 5 J - p ,5r +lir ,+h.c..(10) 

3,(r,r') 

This interaction can be regarded as a part of the 
spin-phonon coupling H sp = -Js P J2j,( r ,r')( S 3,r ' 
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\ J y2(z 2 ) ~ Jy 3 (z 3 )\ and (k x ,k y{z) ) = (tt, tt) for \J Vl ( Zl )\ < 

I ^2/2(^2) — Jys{zs) I- 

We show some typical examples of obtained phase dia- 
grams in Fig. [3j When interchain couplings are not frus- ^ 
trated such as the J yi (J Zl ) dominant case of Fig. [3]^a) and ^ 
(b), the SDW ordered phase is largely enhanced and the 
nematic ordered phase is reduced to a higher-field regime 
compared to the crossover line (K+ = 2) in the J1-J2 
chain. This is because the effective couplings Jfff >W anc ^ 
J^ e are respectively generated from the first- and second- 
order cumulants, and therefore Jffp w is generally larger 
than J^ e in non-frustrated systems in a sufficiently weak 
interchain coupling regime. When both the couplings J y2 



0.004 
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FIG. 3: (color online) Phase diagrams of the 3D magnets ([2]) 
of weakly coupled J1-J2 chains in the M-T plane, which are 
derived by the ICMF approach. The temperatures T S DW(Ne) 
denote the 3D SDW (nematic) transition points. The vertical 
dashed lines denote the crossover lines in the ID J1-J2 chain 
between nematic dominant and SDW dominant TL liquids. 



Sj yr >)(Sj+i^ r - Sj+iy). One can easily expect that 
Eq. (pTU]) enhances the spin nematic ordering. Applying 
the field theory strategy to the system Hsd +%4, we find 
that J^fi is replaced with J^g— 4 J^Co^os k y +cos k z ). We 
thus obtain the phase diagram for Hsd + ^4, as shown 
in Fig. [U Comparing Fig. [3^a) [(b)] and Fig. a) [(b)], 
we see that an inter-chain four-spin interaction definitely 
enhances the 3D nematic phase even if its coupling con- 
stant J4 is small. Since J4 is usually negative, it favors 
ferro-type nematic ordering along the y and z axes, i.e., 



0. 



Conclusion.— We have constructed finite-temperature 
phase diagrams for 3D spatially anisotropic magnets, 
which consist of weakly coupled spin-^ J1-J2 chains, in 
applied magnetic field. Incommensurate SDW and spin 
nematic ordered phases appear at sufficiently low tem- 
peratures, triggered by the critical TL liquid properties 
in the J1-J2 spin chains. We reveal several natures of 
orderings of coupled J1-J2 chains: The 3D nematic or- 
dered phase is generally smaller than ID nematic domi- 
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FIG. 4: (color online) Phase diagrams of the weakly coupled 
J1-J2 spin chains (|2]) with a weak inter-chain four-spin inter- 
action Te- 



nant region, while it can be larger if we somewhat tune 
the inter-chain couplings. The ordering wave numbers 



tend to be or 7r, and a small four-spin interaction 
H4 efficiently helps the 3D nematic ordering. 
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